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Biological macromolecules have complex and non-trivial energy landscapes, endowing them 
a unique conformational adaptability and diversity in function. Hence, understanding the 
processes of elasticity and dissipation at the nanoscale is important to molecular biology 
and also emerging fields such as nanotechnology. Here we analyse single molecule fluctua- 
tions in an atomic force microscope (AFM) experiment using a generic model of biopolymer 
viscoelasticity that importantly includes sources of local 'internal' conformational dissipa- 
tion. Comparing two biopolymers, dextran and cellulose, polysaccharides with and without 
the well-known 'chair-to-boat' transition, reveals a signature of this simple conformational 
change as minima in both the elasticity and internal friction around a characteristic force. 
A calculation of two-state populations dynamics offers a simple explanation in terms of an 
elasticity driven by the entropy, and friction by barrier-controlled hopping, of populations on 
a landscape. The microscopic model, allows quantitative mapping of features of the energy 
landscape, revealing unexpectedly slow dynamics, suggestive of an underlying roughness to 
the free energy. 
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The paradigm underlying many force probe experiments is the Unear increase of tensile force 
on a single biomolecule with time^ An emergent theme from constant loading rate experiments is 
the propensity for conformational change in biomolecules, from reversible processes such as chair- 
boat transitions in polysaccharides^"^ and the overstretching transition in DNA^~^ to the irreversible 
unfolding of concatamers of protein domains^"^"*. In addition, conformational transitions are ubiq- 
uitous in biological processes; for example, static and dynamic changes in structure are known to 
be important in many signalling processes in molecular biology^^. However despite their impor- 
tance, the physical processes that underly these transitions, in particular the role of conformational 
elasticity and internal friction, are poorly understood. 

Despite the success of constant loading rate experiments, they can provide only limited in- 
formation; the elastic response function for each molecule under reversible conditions, and at 
most global dynamical information, such as the rate of unfolding of a protein, from irreversible 
stretching. A case in point is the polysaccharide dextran, which exhibits a reversible plateau in its 
force-extension response, due to a local chair-boat transition that has been shown to be two-state 
in nature^' Such experiments provide the free energy difference and distance between states, 
however, the dynamics of this transition are too fast for stretching experiments to probe. A fuller 
understanding of the response of single biopolymers during forced unfolding or refolding could 
be provided by analysis of the local linear viscoelastic response. Significantly, local dissipation 
would give access to finer-scale conformational dynamics, for example, the rates of transitions be- 
tween different states along the unfolding or refolding pathways of a protein. A close analogy is 
found in the macroscopic rheology of complex fluids, whose dissipative mechanical spectra reflect 
dynamics of various structural, molecular and topological transitions^'. 

Recent experiments^^"^^ measuring the viscoelastic properties of single biomolecules as a 
function of force, including polysaccharides and proteins, have gone some way to achieving this 
goal. The results show highly non-trivial features, particularly in the dissipative part of the spec- 
tra, where measured frictions are many orders of magnitude larger than that due to solvent. In 
particular, the friction of dextran exhibits a minimum at a force that coincides with the plateau in 
the force-extension trace, indicating it arises through a process related to the local conformational 
transitions in the chain^^'^"^. In addition, although it is clear that a plateau in the force-extension 
response, should give rise to a minimum in elasticity, the underlying statistical mechanics of this 
change are not well understood^' ^. Here we seek to understand the origins of these features in the 
viscoelasticity of dextran and by doing so give broad insight to the nature of elasticity and friction 
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for simple conformational transitions. 

Dextran and cellulose are polysaccharides that are biological polymers composed of glu- 
cose monomers, a six-membered ring molecule, which is known to have a number of stable 
conformations^^ (Fig.l). These biopolymers differ by the way the glucose ring is linked into the 
backbone of the polymer. In dextran, one of the linkages is axial to the plane of the ring and 
thus force promotes conversion from the nominally stable chair state to a more elongated boat- 
like conformation, in which this linkage is equatorial^^, as shown in Fig. la. This gives rise to 
dextran's characteristic plateau in its force extension response (Appendix: Fig.l). In contrast, the 
glucose ring in cellulose is already near maximum elongation since all its linkages are equatorial 
to the plane of the ring (Fig. lb) and results in almost ideal freely jointed chain (FJC) properties 
in its force-extension response^' (Appendix: Fig.7). Hence, we will show through the experi- 
mental comparison of the viscoelasticity of these two polysaccharides, that the two-state nature of 
the transition in dextran and its absence in cellulose, provides an ideal test-bed to understand the 
characteristic viscoelastic response of simple forced conformational transitions. 

We determine the viscoelasticity of dextran and cellulose, using a recently developed tech- 
nique for measuring the Brownian dynamics of single molecules under force-clamp conditions^. 
Fig.2a summarises the experimental apparatus and procedure, with details given in Methods. The 
principle of the experiment is to hold a single molecule between tip and substrate of an AFM at 
constant force, whilst observing the thermal fluctuations of the cantilever. The fluctuations con- 
tain inherent viscoelastic information, which we obtain via calculation of their frequency power 
spectral density (PSD). A conventional proportional-integral-derivative (PID) feedback loop with 
a response time of ~10 ms, monitors the cantilever deflection signal and adjusts the piezo substrate 
to maintain a constant average force (F) or "force-clamp" on the molecule between the tip and sub- 
strate. A key idea of this technique is that by controlling the force we probe the local viscoelasticity 
of single biomolecules as they explore their energy landscape under near equilibrium conditions. 
Measurement of the force-dependent power spectra is exemplified in Fig. 3 for cellulose, where 
it is clear that the PSD peak position, width and amplitude are dependent on the response of the 
biopolymer. 

To quantify these changes and extract viscoelastic information from the thermal spectroscopy 
power spectra, we model the biopolymer using a modified Rouse model that includes local con- 
formational internal friction, in addition to solvent friction^^'^^. The Rouse model is a generic and 
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highly successful description of the coarse-grained dynamical behaviour of polymers^"' where 
we note that in the typically highly extended conformations in our experiment, long-range hy- 
drodynamics give only logarithmic corrections to local drag. The Rouse with Internal Friction 
(RIF) polymer is represented as a series of beads with solvent friction (so, connected by spring and 
dashpots of elasticity kq and internal friction Qq. In the continuum limit, internal friction adds an 
extra term in the standard Rouse equation, which describes a dissipative force proportional to the 
rate of change of local conformation, represented as the coarse-grained curvature of the chain. 



^ dR{n,t) f ^ d\d^R{n,t) , 



where R{n, t) represents the space curve of the polymer with contour variable n, subject to a 
local Langevin force f(n,t), which is uncorrelated for different times. Normal mode solutions 
of this equation decay in a single exponential manner with a mode dependent relaxation time, 
Tp = 27r^fop2 , where p is the mode number, the effective mode friction is renormalised compared 
to standard Rouse theory to (p = 2N(^sq + ^ , where intuitively, the new term accounts for 
an increasing internal friction of higher curvature short wavelength modes. AFM experiments 
probe the end-to-end vector of the polymer, whose response can be found by summing over all odd 
modes; in frequency space this gives the following useful closed form expression for the dynamic 
compliance: 



tanh ? 



2N V ^ V 1+*'^'^* 
Jar{^) = j====^, (2) 



where tr = A^^Cso/tt^ko is the contribution to the relaxation time of the first mode due to solvent 
friction and Tj = / f^o is the mode-independent contribution to the relaxation time due to internal 
friction. This model successfully encompasses the behaviour of both types of friction; in the limit 
of large internal friction (tr « Ti), Eq.© reduces to a single mode spring and dashpot model, 
JAni^) = (i+iW ) ^^'^ when solvent friction dominates (tr >> r^) to the Rouse model, given 
by the limiting form, Jar{^^) = (|-\/za;r/j) / y/itoTR, till a critical frequency l/xj, when 

the internal friction of high curvature modes dominates to give single mode relaxation again. The 
Fluctuation-Dissipation Theorem (FDT)^^, P{lu) = —2kBTJ"{uj)/uj is then used to calculate the 
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total power spectrum P(a;) of a RIF polymer combined with a SHO response of the cantilever and 
cantilever, where J"{oj) is the imaginary part of the response function J{uj). 

Shown in Fig.4a are the effective monomer elasticity of cellulose and dextran, from the RIF 
model fits and normalised by contour length. In previous work^"* we showed that calculating the 
elasticity spectrum directly from the numerical derivative of extensible FJC fits, and secondly from 
the thermal spectroscopy method agree very well. We verify that using the more refined RIF + 
cantilever model to analyse the PSD, also provides very good agreement. As previous studies have 
shown^"^'^^'^^'^'^^, at low force (in these experiments), elasticity is due to the reduction of chain 
conformational entropy as it approaches its contour length, after which contour length elongation 
with constant elasticity becomes more favourable. At higher force, however, the minimum in the 
elasticity spectra for dextran at ~ lOOOpN, which is absent in the cellulose spectrum, marks a clear 
signal of the conformational transition in the former. 

The key advance afforded by using the RIF model in analysing the PSD is the new informa- 
tion about the two sources of dissipation, not distinguished in previous work^^'^^'^"*'^^; the solvent 
friction and internal friction of the single biomolecule. We find consistently from the RIF analy- 
sis, that solvent friction is very small within the errors of this experiment (< 0.01//gkHz). Hence, 
these chains are 'short', as defined by N -x/Cio/Cso where N 400, which indicates that 
dissipation is dominated by internal friction at high stretch, and explains the success of the spring 
and dashpot model in previous modelling of the dissipation of dextran 18,19,24,25 -pj^g internal 
friction force spectrum itself exhibits non-trivial behaviour as shown by the comparison of cellu- 
lose and dextran in Fig.4b. At low force, both polymers show an increasing internal friction with 
force, followed by a plateau. Crucially, at higher forces, the spectra of cellulose and dextran dif- 
fer; qualitatively, the minimum in the internal friction force spectrum of dextran at ~ lOOOpN and 
its absence in cellulose, confirms that source of this change in the friction of dextran is from the 
chair-boat conformational transition of the glucose ring. 

To make this conclusion more concrete we link the features of the experimental elasticity 
and friction force spectra to the conformational transition in dextran, using a simple model of 
population dynamics on a discrete 2-state energy landscape, which we show predicts the same 
viscoelastic signature of simple forced transitions, as seen in Fig.4. The parameters of the dis- 
crete 2-state model are as described in Fig.Sa, in which we assume populations obeys Boltzmann 
statistics and dynamics follow activated Arhennius transition rates. Using an approach similar 
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to ' , the effective response of the populations at a frequency uj can be calculated by applying 
an oscillatory force /o cos ujt to the energy landscape. The dynamics of the population p{t) in 
state 1 (say, the short state, so that probability of extended state is 1 — p{t)) are then described by 
^ = — (Ai2(t) + \2i{t))p{t) + A2i(t), where the rates A12, A21 vary with time due to the oscillating 
perturbation of the landscape. In the Brownian linear response limit {f^x <^ ksT), there are in- 
phase and out-of-phase oscillating solutions to this differential equation, such that the extensional 
response of the monomer is a simple spring and dashpot Jui^) = ki2+<i2<^ ' which we identify 
the effective elasticity and friction as: 
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Ci2(F) = ^(ri2(F) + r2i(F)), (4) 

and where Pq{F) = (1 + e"^^*^*^^^)"^ is the equilibrium Boltzmann probability for the short state. 
In addition, T12 is the forward hopping time between states and T21 is the corresponding backward 
time, where in general r^j = tqc^^'^'j'''^^ with j3 = l/ksT, AGJ^ the free energy barriers for 
interconversion and tq = 27r(b/ f^b is a prefactor that arises from mapping the Kramers' first passage 
problem on a continuous free energy landscape^^ G{x) to a discrete description (Fig. 5a), where Cb 
and Kb are the effective friction and curvature of the barrier. 

Plotting these (Fig.5b&c - on a natural logarithmic scale to emphasise their exponential na- 
ture) we see a characteristic minimum in both the elasticity and internal friction force spectra. In 
the former case, it is clear that the source of the elasticity is entropic in nature and not enthalpic 
as has been previously asserted^; force controls the shape of the energy landscape or the relative 
populations of monomers in short or extended states and hence, the effective 'size of box' that the 
monomer can explore. So Eq.© is an expression of the equipartition theorem k = ksT / (A6^), 
where (A6^) is the mean square fluctuations of the monomer; in Fig. 5b at low force, AG{F) is 
large and positive, hence monomers are confined to the short state, fluctuations (A6^) are small 
and the effective stiffness is large. As force decreases the energy difference, populations spread 
across the two states and the effective size of the box (A6^) increases, causing the stiffness to de- 
crease (exponentially). The stiffness subsequently passes through a minimum when AG{F) = 
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and (A6^) is maximum, corresponding to a state of maximum entropy, when the probabilities to 
be in either of the states are equal. On further increase of force, AG{F) becomes negative and 
monomers become increasingly confined to the extended state ((A6^) decreasing) and the stiffness 
increases exponentially. It is simple to see that the elasticity is purely entropic, since any enthalpic 
contributions to the free energy difference AGq can contribute only linearly to the free energy as 
the extension of the monomer is increased. Thus, the molecular elasticity of a monomer is defined 
by its entropy on a discrete energy landscape. 

Eq.® predicts that the internal friction for a 2-state landscape is proportional to the sum of 
the times to interconvert from state 1 to 2 and back, from state 2 to 1 . Applying a force to the 
monomers changes the activation barriers to interconversion, which changes the average time to 
interconvert and thus ultimately, the internal friction. Fig. 5c shows schematically how the inter- 
nal friction should vary with force in a discrete two-state landscape. As force lowers the barrier 
AG\2{F) of interconverting from 1 — 2, the internal friction should decrease, passing through a 
minimum when the barriers on either side are approximately equal (when xi = X2 this occurs at 
exactly AGfa = AG^i) and then increase again at high force as the barrier for the reverse transi- 
tion (AG|]^) and hence T21, becomes large. It is interesting to note that, whilst the hopping time 
passes through a minimum, the corresponding relaxation time r = r^^^ + must pass through 
a maximum, since relaxation is dominated by the smallest barrier. Hence, on average fluctuations 
away from equilibrium occur on a hopping timescale r* ~ T12 + r2i, whilst relaxation back to 
equilibrium occurs on the timescale r. We see Eq.©, is a microscopic fluctuation-dissipation re- 
lation for a discrete bistable landscape, that links friction to the timescale for fluctuations due to 
activated barrier-hopping. 

Useful information about the position of the transition state can also be obtained by analysing 
the relative positions of the minima in the elasticity and internal friction spectra. The difference in 
the forces at which the minima occur AF, can be found from the derivatives of Eq.@ and Eq.©: 



and thus provides information on the relative position of the transition state, xi or X2 (Ax = 

Xi + X2). 




(5) 
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The results of this population dynamics model, thus provide a simple way to understand 
the minima in the elasticity and internal friction spectra, in terms of entropy and barrier- hopping. 
However, to understand the entire force regime (~ 100 ^ ISOOpN), in addition to the viscoelas- 
ticity of the 2-state conformational transition, we need to include the physics of the chain at low 
and intermediate forces, before the critical force at which the conformational transition occurs. At 
low force we use a Frictional Freely Jointed Chain (FFJC) model (Appendix) of rods intercon- 
nected with joints with constant friction Qe to give an elasticity Kpjc{F) = and an internal 
friction that increases linearly with force, Cfjc{F) = 2*^-^' which are both valid at high stretch 
(F ^ ksT/b ~ 4pN for b ~ Inm). At intermediate force we account for the local viscoelastic- 
ity of stretching a dextran monomer in the short or extended states, using constant elasticities ki, 
K2, and internal frictions (^i, C2, respectively. We assume that these processes add mechanically 
in series, since they provide independent and additive extensions to the overall chain length (see 
Methods). 

Fitting to the elasticity force spectra of cellulose and dextran (normalised by contour length), 
we find excellent agreement as shown in Fig. 4a, where the solid line represents the full elasticity 
Eq.© generated using the average of the parameters determined over a number of single molecule 
experiments (Cellulose: ki = 36000 ± 18000 pN/nm, b = 1± 0.5nm; Dextran: AGq = 16.5 ± 
OAksT, Ax = 0.066 ± O.OOSnm, ki = 10000 ± lOOOpN/nm, K2 = 39000 ± 2000pN/nm, b = 
0.63 ± 0.02nm). These values agree well with the literature 2.3,5,38 confidence we can 

describe the whole elastic force spectra for both cellulose and dextran; at low force (below 800pN) 
stiffness increases as entropy is lost due to the orientation of monomers along the line of force and 
finally reaches a plateau representing a constant stiffness due to the enthalpy of stretching the bonds 
comprising the glucose ring. However, the response of dextran differs dramatically at higher force 
as the more extended state becomes thermodynamically favourable. Within the framework of the 
2-state model presented, the subsequent decrease in stiffness can be understood since it becomes 
more entropically favourable for the chain to elongate. Interestingly, in the elasticity spectrum 
of dextran, at around 400-500pN, the model slightly, but consistently, underpredicts the elasticity 
below this force and overpredicts it above this force. This plateau may be explained by the entropic 
elasticity of other internal states, possibly the C5-C6 bond rotation in dextran 

In performing fits to the internal friction spectra, all elastic parameters are constrained to val- 
ues obtained from fits to the elasticity spectra (see Methods). Below we discuss quantitative values 
of each of these friction processes separately, even though actual fits were performed globally 
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across the whole force range. 

Firstly, we examine the effective internal friction associated with stretching the glucose 
monomers in their various conformations. For cellulose, we find (^i = 110 ± 50/igkHz and for 
dextran, = 25 ± lO/xgkHz and ^2 = 120 ± 50;UgkHz, for the short and extended states, re- 
spectively. Strikingly, these numbers are roughly 7 orders of magnitude larger than the friction 
expected due to solvent {( = Qixrih ~ 10~^/igkHz for h ~ Inm). The most plausible source for 
such a high local effective friction is roughness in the free energy landscape. A model of dynamics 
on a rough Gaussian landscape with RMS energy fluctuations e:"^" predicts a sensitive enhancement 
to the effective friction constant C* = C 6xp {e/ksTY, giving an effective roughness for stretching 
cellulose and dextran as £ ~ AksT. For comparison, recent constant loading rate experiments'^^ on 
the protein imp-/3, using theoretical results in"^^ suggest a Gaussian roughness of order e ^ S.T/cbT. 
In the case of these polysaccharides, this roughness may arise from the many sub-states separated 
by barriers that must be traversed in stretching the monomers; for example; for example, there 
are many conformations of a glucose ring, (in total 14 canonical chair, boat and twist-boat confor- 
mations, separated by 12 half-chair and 12 envelope conformational transition states^^ "^^), which 
will contribute to extension and may become more or less favourable under tension. In addition, 
the hydroxyl groups of glucose give rise to the possibility of intra and intermonomer hydrogen 
bonding, as well differing degrees of solvent accessibility. Such states are particularly suggested 
by slow undulations in the elasticity and internal friction spectra of cellulose for forces greater than 
1000 pN. 

At low force, we see a similar picture for the 'joint' friction of the FFJC model, obtaining 
values of the order (g ~ l^ugnm^kHz (cellulose: (g = 0.9 ± 0.7/ignm^kHz; whilst for dextran 
errors from fits suggest Q < 1.2yUgnm^kHz). These numbers are roughly 6 orders magnitude 
greater than the friction of a rod of length b rotating in a solvent {7rr]b^ /A ~ 10^^/ignm^kHz). We 
can again appeal to an underlying molecular explanation, where joint friction is due to hopping 
between dihedral angular states, with an average hopping time of Thop ~ Ce/knT ~ 0.25msec. 
Again, these very slow dynamics are suggestive of an underlying roughness to the rotational free 
energy (e ^ 3.7A;bT, where C7C ~ 10^). 

In the case of dextran, the marked decrease in internal friction around ~ lOOOpN, contains 
information on the dynamics of interconversion between the short and extended states, for which 
Eq.© provides a simple model. In principle, fitting to the internal friction spectra would determine 
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the position of the barrier xi (with constraint X2 = Ax — Xi) and the zero-force interconversion 
times Ti2{F = 0) and r2i(F = 0) (measurements at different temperatures could in principle, 
determine the free energy barriers for conversion in each direction). However, although we find 
very good fits around the transition region, they are underdetermined, due to the low frequency 
restriction of the data. 

To constrain our fits further, we use Eq.©. Inspecting Fig.4a&b (solid squares), indicates 
that AF ^ ± lOOpN, given that the spacing of points in the spectra is approximately lOOpN. 
However, negative values of AF imply from Eq.© a transition state that is closer to the short state 
than long (— < 1), which is not feasible on geometric grounds, given that its curvature is roughly 
four times smaller than the extended state (^ ~ |) and that the forward free energy barrier must 
obey AGq^2 > AG'o(= IQ.bksT). Hence, a reasonable assumption is that < AF < lOOpN, 
implying a position of the transition state in the region 0.033 < Xi < 0.053nm. Fitting to the 
friction spectra, so as to satisfy this constraint on xi, we find lns< r2i(0) < 100ns and 0.01 
s< ri2(0) < 1 s (see Methods). From Kramers' theory of activated diffusive barrier crossing, 
the exponential prefactors of these times are related to the curvature and friction (b of the barrier, 
which when mapped onto a discrete landscape is given by tq = 2TxCb/ i^b- Thus, given an order of 
magnitude estimate of the barrier friction C,}, ~ |(Ci + C2) ~ 70/igkHz and that tq < r2i(0), 
we find that the barrier must be very sharp; given by the following approximate bound, > 
lO^pN/nm. Fig. 6 shows a graphical to-scale reconstruction of the free energy landscape based on 
the parameters extracted from the modelling of the viscoelastic force spectra of dextran. 

In summary, we have shown how macroscopic ideas of elasticity and friction can be extended 
to the study of the energy landscape of conformational transitions. Eq|3]and Eq|4|are in essence 
microscopic equivalents of the equipartition theorem and the diffusive fluctuation-dissipation re- 
lation, where the spatial and temporal properties of the fluctuations are determined by the shape 
of the energy landscape, which in turn determine its effective elasticity and friction. In the case 
of dextran, applying tension to its energy landscape, drives an entropic transition, where elasticity 
and friction decreases as the populations become more spread and barriers are lowered. These 
ideas are of wide relevance; from applications to the field of molecular nanotechnology, where 
microscopic processes of elasticity and internal friction may guide and constrain engineering de- 
sign, to understanding fundamental processes of molecular biology, such as the study of internal 
transitions in biomolecules, including the action of molecular motors, allosteric signalling, force- 
sensing between cells, stretching transitions in DNA and RNA, and emerging data on elasticity and 
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dissipation from the fluctuations of a refolding protein. 
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Methods 

Experimental Materials & Methods The protocol used for thermal force-clamp spectroscopy 
is as described in^'*, we summarise the procedure here. The first part of the experiment follows 
conventional force-spectroscopy protocol, where the cantilever is pressed into a polysaccharide 
monolayer with a force ^ 10 nN for ~ 1 s, after which it is retracted from the substrate at a 
constant speed. When a pre-determined force set-point is reached, the force-clamp protocol is 
initiated, which involves either reducing force in discrete steps of ~ 100 pN and being held for 
~ 3 s, or reducing force slowly and continuously at ~ 8 pN/s. In some measurements we have 
used this latter continuous approach, however, both procedures produce the same results within the 
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errors of each method (not shown). In either method the force is controlled using a proportional- 
integral-derivative (PID) feedback loop with a response time of ~ 10 ms, whereby the cantilever 
substrate separation is adjusted to maintain a certain cantilever deflection. A response time of 
~ 10 ms means the feedback loop cannot respond to fluctuations faster than 10 ms. Thus, for 
frequencies greater than ~ 0.1 kHz, an average force F is maintained. After the force-clamp phase, 
the cantilever is again retracted from the substrate at a constant speed, till at some critical force 
the polymer detaches. Immediately after detachment, the PSD of the free cantilever is recorded 
as the cantilever is brought towards the substrate in 30 nm steps. These free cantilever PSD are 
then fit using a simple harmonic oscillator model (SHO), Pc{io) = -m!^'^)^\cc'^'^ ' '^^^^^^^^S '^he 
cantilever effective stiffness friction constant Cc and mass rric. These parameters then serve as 
constraints in the curve fits to the power spectra of the cantilever/molecule system. 

To extract the elasticity, internal and solvent friction as functions of force, we treat the force 
clamp experiment as two linear system elements in parallel, since the change in extension of the 
polymer and cantilever are the same at their point of contact. It can be shown (Appendix), that for 
a system in parallel the total dynamic compliance of the system Jt{^^) is given by 



where Jx {oj) is the dynamic compliance of the cantilever, for which we use a SHO model {Jx (i^) = 
{k — muj'^ + iCuj)^^). This is just the frequency-dependent extension of the parallel addition 
that arises naturally in our experiment. We then use the fluctuation dissipation theorem (FDT) 
P{ijj) = —2kBTJ"{uj)/uj (J" represents the imaginary part of the complex function J) and 
Eq.(0) & Eq.® to calculate the total power spectrum of the cantilever -i- RIF polymer. This enables 
measurement of the elasticity, internal friction and solvent friction as functions of force, for exam- 
ple, as shown by the RIF model fits to the power spectra of cellulose in Fig. 3. In fitting the PSDs to 
the RIF-i-cantilever model, we constrain the chain solvent friction to be between < < 0.1 
/igkHz, since a reasonable estimate of the solvent friction is given by QixriLcapproxA x 10"'^ 
/igkHz, where the contour length Lc ~ 200 nm, is the typical contour length of molecules in 
the experiment, representing the maximum effective hydrodynamic radius of the chain^"^. 

Fitting to Elasticity and Internal Friction. In both the elasticity and internal friction force spec- 
tra, there are a number of different physical processes that underly the observed behaviour. A 
reasonable assumption is that the noise on each physical processes is uncorrelated, so that the total 
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power spectrum is the sum of the power spectra of each process. In the low frequency regime of the 
experiments (lot <S 1), the rules for summing the elasticity and frictions of the different processes 
are then: 

<F) = (K^'jciF) + Ki2iF)+Po{F)K^' + {1-Po{F))k^Y' , (7) 

where importantly, the bond elasticities of the short and extended states are weighted by the prob- 
ability to be in those states at a given force, where Po{F) = (1 + e^^'^(^))^. In fitting to the 
internal friction force spectra, we use the parameters extracted from fitting to the elasticity spectra 
as a constraint to the fits. Through trial and error with fits with different fixed values of r2i(0), we 
found the values of the zero-force backward hopping time that correspond to the bound calculated 
on xi in the main text. We checked that the values of Ti2(0) also determined from the fits, were 
consistent with detailed balance (ti2(0)/t2i(0) = 6^^^^°). 
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Figure. 1. Structure of dextran and cellulose, (a) Simplified diagram of the molecular structure 
of dextran, which is an a-(l— i>6) linked polysaccharide of glucose, where the monomer length is 
defined by the distance between adjacent non-ring oxygens on the backbone, as shown schemati- 
cally. The a linkage at CI is axial in the lowest energy ^Ci chair conformation^^ (above), which 
under application of a tensile force-field promotes one of a number of more elongated boat or 
skew-boat conformations, of which the ^'^B is shown (below)^'^^ "^^ '^'^. The increased length Ax 
gives rise to a plateau in force-extension measurements (Appendix Figure 7) as the more elongated 
boat-like conformations are populated under increasing force, (b) Cellulose on the other hand is 
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a /?-(!— i>4) linked polysaccharide of glucose, whose equatorial linkage at CI in the chair state, 
means the monomer is already near maximum elongation and a its force-extension behaviour fol- 
lows simple polymer elasticity models due to reduction of chain entropy at high stretch (Appendix 
Figure 7). 

Figure. 2. Force-clamp thermal noise spectroscopy, (a) Schematic diagram of the experimental 
setup for thermal noise spectroscopy. (b)&(c) show the force clamp protocol used: (b) typical 
experimental force-extension traces showing a retract, approach and retract cycle, using dextran 
as the sample polymer, where the traces have been offset for clarity. The characteristic shape of 
the final curve (curve 3), which in the case of dextran exhibits a shoulder indicative of the well- 
known conformational transition in dextran^, confirms that only a single molecule was attached, 
(c) experimental force-time trace, where the numbers and colours correspond to the same sequence 
in (b). The force-clamp phase (phase 2) lasts for a total of 42 seconds, where the dextran polymer 
is held for 3 seconds at each of 14 discrete forces (the initial and final extensions, phases 1 & 3, 
are shown on an expanded time scale). 

Figure 3. Force dependent PSD of single molecule of cellulose. Comparison of the PSD of fluc- 
tuations of cantilever tip, when free (black circles) and with a single cellulose molecule attached, 
held with forces of 320 pN (green), 620 pN (purple) and 920 pN (red). The solid lines correspond 
to fits using either a simple harmonic oscillator model for the cantilever (black), or the RIF model 
of the biopolymer combined with the cantilever (green, purple and red solid lines) described by 
Eq.© & Eq.©. 

Figure 4. Viscoelastic force spectrum of cellulose and dextran. (a) Elasticity force spectrum and 
(b) internal friction force spectrum multiplied by contour length of each molecule (giving the 
inverses of the compliance and mobility per unit length) for cellulose (solid diamonds) and dextran 
(solid squares, where Lc is obtained from FJC fits to their respective force-extension traces (Ap- 
pendix Figure 7)). Data points represent measurements using thermal force clamp spectroscopy, 
where different colours represent separate single molecules. The solid lines represent curves gen- 
erated using the full elasticity (a) and internal friction (b) expressions given in Eq.© & Eq.®, 
using the average of the parameters determined over all of single molecule experiments (see main 
text), apart from Qe = 0.6yUgnm^kHz (half the upper bound in the main text) and 2-state internal 
friction parameters derived in the text consistent with a fixed zero-force backward interconversion 
time r2i(0) = 100ns (i.e. Axi = 0.053nm, ri2(0) = Is). Horizontal error bars represent an ap- 
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proximate 10% systematic error between experiments in determining the true force scale, through 
errors in measuring cantilever elasticity and cantilever deflection sensitivity. Vertical error bars 
represent errors from the fits to the PSDs. 

Figure 5. Viscoelastic Force Spectrum on a discrete bistable landscape, (a) Schematic di- 
agram of the discrete free energy landscape used to calculate the elasticity and internal fric- 
tion force spectra (Eq.@ & Eq.©). (b) Elasticity force spectrum on a discrete 2-state land- 
scape; force controls free energy difference AG{F) = AGq — FAxi2 and hence spread (A6^) 
and elasticity ki2{F) = ksT / {MP')). Elasticity is entropic in nature as elasticity decreases 
in direction of increasing entropy of monomers, (c) Internal friction force spectrum for a dis- 
crete 2-state landscape; force controls activation barrier heights {AG\2{F) = AG1^2 ~ FAxi, 
AG2i(-^) = ^^021 + ^^^2), and therefore also the internal friction. Hence, at a given force, 
internal friction is dominated by the activation barrier that is largest. 

Figure 6. To-scale reconstruction of continuous free energy landscape of glucose, based on pa- 
rameters extracted from theoretical modelling of the viscoelastic force spectra of dextran. Dashed 
features indicate areas of landscape that are uncertain, for example position of barrier, or infor- 
mation unattainable with current experiments like the activation barrier heights. Barrier curva- 
ture shown is K5 = lO^pN/nm. Grey lines indicate a roughness to the landscape with RMS 
deviation e ~ AksT, as a plausible interpretation for significantly enhanced friction of wells. 
AG[) = AGo + kBT\n{^Ki/K2) = (16.5 - \n2)kBT ^ IS.S/cbT is the free energy difference 
between the minima of a continuous landscape, which excludes the entropy of vibrations of the 
wells. 
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Appendices 

Dynamic compliance of parallel system elements Here we derive the total dynamic compliance 
or response function Jt{oj) for the cantilever and polymer in parallel, which each have response 
functions Jx((^) and JAii(<^), respectively. Starting in the time domain, we can write down the 
solution for the cantilever and polymer motion as: 



where F represents an external force applied to the system and / the internal force that they share 
according to Newton's 3'"'^ Law. By definition, the Green's response of the whole system is its 
response to a unit impulse of force, so we let F{t) — rjS (t), where 6 (t) is the Dirac "delta-function" 
and T] the size of impulse. This gives 
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AR{t) = lj^R{t) - f J^nit - t')f{t')dt' 
^ Jo 

x{t) = px{t) + r Jx{t - t')f{t')dt' 

^ Jo 



Taking the Fourier Transform of these (presuming all response functions are zero for t < 0) we 
find 



ARico) = I Ja/?(^) - JAR{uj)f{uj) (1) 



X(o;) = |jxH + JxM/M (2) 

Thus, using the fact that the cantilever and polymer displacements must be the same for all times 
(Ai?(t) = X(t)), we can solve for the internal force f{uj): 



2 JarH + Jx{to) 



this can then be plugged back into Eq.([T]) or Eq.©, to give the total dynamic compliance as the 
displacement response due to a unit delta function input: 



^ , , AR{uj) X{uj) Jx{uj)Jar{uj) 

Jt\^) — 



rj T] Jx{uj) + Jar{uj) 



Frictional Freely Jointed Chain To model the molecular viscoelasticity of a polymer at forces 
which are small (approximately, F < 500pN), we develop a Frictional FJC (FFJC) model of rods 
interconnected with joints with constant friction and calculate the form of (pjc{F). We focus on a 
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single monomer, assuming that each rod of the FJC is statistically independent, so the stiffness of 
each rod will add mechanically in series to the stiffness of the whole chain. Typical monomer/rod 
lengths for polysaccharides are b ~ Inm, so our experiments are in the regime where F ^ ksT /h 
and the elasticity spectrum can be calculated from statistical mechanics as 



.„c(F) = ^ (3) 

To model the internal friction of a FJC we again focus on a single monomer/rod in the high force 
regime and consider that to rotate such a rod there is some friction C0 opposing this motion, which 
we presume is constant and associated with the internal friction of 'joints' between rods. The 
rotational equation of motion for a segment or rod of length h held under a large tensile force F 
(Fb > kBT) will be 



(em = -FhO (4) 

Now we consider how these dynamics project onto the line of applied force. The change in pro- 
jected length of the monomer compared to its actual length will be A6 = 6(1 — cos^^), which in 
the small angle limit will be: 



A6 ^ -6^2 (5) 
2 ^ ' 



Differentiating A6, and using Eq.© & Eq.©, we find its equation of motion to be: 



Ab = bee = — = — —Ab 

Ce Ce 



This again has an exponentially decaying solution Ab{t) ~ e ^/'^fjc^ with time constant Tpjc 

Ce 
2Fb 



Now Tpjc = Cpjc/f^Fjc and thus, using Eq.©, the effective friction along the z direction is 
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then 



which predicts a linear increase of the internal friction constant with force, presuming (0 is con- 
stant. 
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Fig. 7. Extension-Force traces from constant pulling rate experiments for (A) 4 different dextran 
molecules and (B) 5 different cellulose molecules. Each trace represents the final stage (marked 
(3) in Fig. 2 in main text) of the thermal noise force-clamp spectroscopy protocol from which we 
determine the contour length Lc. For cellulose we fit the extension as function of force, using a 
FJC of elastic segments (FJC-i-), where the chain is characterised by a number of Kuhn segments 
Nk, which have length b and elasticity n to give a contour length Lc = Nkb. For dextran, we 
assume a Boltzmann weighted sum of the two states, where the monomer lengths in the chair and 
boat state are represented by different Kuhn segment lengths and elasticity of a FJC-i- model: 

where C is the Langevin function, F is the imposed tension, and AG{F) = AGq — FAxi2, with 
Axi2 being the spatial separation between the two states. Factors of fc^T are dropped for clarity. 
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For dextran, we assume the contour length is given by ~ Nj^hi. 
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